Geometry and topology of singularities in spherical dust collapse 
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We derive some more results on the nature of the singularities arising in the collapse of inhomo- 
geneous dust spheres, (i) It is shown that there are future-pointing radial and non-radial time-like 
geodesies emerging from the singularity if and only if there are future-pointing radial null geodesies 
emerging from the singularity, (ii) Limits of various space-time invariants and other useful quantities 
(relating to Thome's point-cigar-barrel-pancake classification and to isotropy/entropy measures) are 
studied in the approach to the singularity, (iii) The topology of the singularity is studied from the 
point of view of ideal boundary structure. In each case, the different nature of the visible and 
censored region of the singularity is emphasized. 
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I. INTRODUCTION AND SUMMARY 

The role played by regular initial data in determining 
the final state of marginally bound spherically symmet- 
ric collapse of inhomogeneous dust is well understood. A 
comprehensive description of this role was first given in 
jl| where explicit necessary and sufficient conditions on 
the initial data were derived for the existence of a naked 
singularity. A naked singularity here means a singular re- 
gion from which there emanates a future-pointing (f.p.) 
causal geodesic. The result referred to dealt with the 
existence of f.p. radial null geodesies. Motivated by a de- 
sire to (i) study the stability of this result with respect 
to the introduction of angular momentum and (ii) gen- 
eralise, in the relevant subsets of the initial data space, 
to all f.p. causal geodesies, we have studied non-radial 
null geodesies in spherical dust collapse 01. We showed 
that there exist f.p. non-radial null geodesies emanating 
from the singularity if and only if there exist f.p. radial 
null geodesies emanating from the singularity. It would 
be useful to know whether this result can be extended to 
timelike geodesies in order to give a complete account of 
all possible causal geodesies. Interesting results in this 
direction were recently given in [0. 

The existence of non-radial geodesies outgoing from 
the singularity might suggest that, topologically, the sin- 
gularity is not a point. It is therefore interesting to study 
in more detail the structure of the singularity and in par- 
ticular to investigate its topological properties. This is a 
difficult problem since, firstly, the singularity is not part 
of the space-time and there is no unique definition for 
the space-time boundary (see for a review). Secondly, 
given the complexity of the geodesic equations, it is diffi- 
cult to study the topology of a given singular boundary. 

In order to have a first intuition about the singularity 
'shape', we can start by studying the behaviour of the 
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expansion of dust world lines in the neighbourhood of the 
singularity. In this way, we can classify the singularities 
as points, cigars, barrels or pancakes ||. 

However, this classification says little about topolog- 
ical issues. For example, in a dust filled Friedmann- 
Lemaitre- Robertson- Walker universe the big bang singu- 
larity is point-like in this classification, but is foliated by 
2-spheres according to a scheme used below. Further- 
more, there are simple examples where this classification 
scheme does not give any results such as for the singu- 
larity in Schwarzschild space-time with negative mass. 
We recall that in the context of the self-similar collapse 
of a scalar field, Christodoulou || has found an exam- 
ple a singular boundary (see below) which is foliated by 
points in some sectors and by spheres in other sectors. 
So, a space-time singularity might not have a distinct 
unique topology and can possess complicated geometri- 
cal features (see also e.g. the Curzon singularity 0). It is 
therefore of importance to investigate some of these non- 
trivial issues in the context of spherical dust collapse. 

This paper aims to study carefully some aspects of the 
singularity and its topological properties in Lemaitre- 
Tolman (LT) space-times. In section II we briefly de- 
scribe the marginally bound LT model and recall previ- 
ous results of pj and H. In section III, we consider the 
problem of extending the result of 0J to all f.p. causal 
geodesies. In section IV, we investigate the causal struc- 
ture of the singularity by studying the behaviour of in- 
going radial null geodesies. In section V, we shall use 
Thome's classification scheme, which turns out to be of 
interest in this case beyond merely classifying the singu- 
larity. We shall also study the behaviour of a number 
of space-time invariants in order to investigate the struc- 
ture of the singularity. These invariants involve the shear 
and the space-time curvature and have been used e.g. in 
studies of the Weyl curvature hypothesis |8| as well as in 
the context of isotropic singularities ||. In section VI, we 
address topological issues by considering non-radial null 
geodesies emerging from the singularity for self-similar 
solutions. We end with some conclusions and comments. 
All our considerations are restricted to the local visibil- 
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ity of the singularity. Another important point to keep in 
mind is that when dealing with radial geodesies, we will 
for the most part treat them as curves in the Lorentzian 
2-space ^constant, <p ^constant. So a 'unique' radial 
geodesic as referred to in e.g. Proposition 4 is a unique 
curve in the r — t plane, but is a 2-parameter family of 
curves in the space-time. The exception is Section VI, 
where the different angular parameters of radial curves 
play an essential role. A bullet • indicates the end of a 
proof. We use 8irG = c = 1 . 



that only the subset {(r, t) — (0, tp)} (where to := t c (0)) 
of this singularity may be visible [pL4|| . Assuming a dust 
sphere of finite radius, we can restrict the range of r to 
[0, b] for some b > 0, and match it to a Schwarzschild ex- 
terior. During the collapse of the dust sphere there can 
be also a curvature singularity given by R r — along 
t = t sc (r), where 



t sc (r) = 2- 



(0) 



II. FIELD EQUATIONS AND SINGULARITIES 

We study spherical inhomogeneous dust collapse for 
the marginally bound case. The line element is [10[ [ll], 
H2l 



ds 2 = -dt 2 + (R r dr) 2 + R 2 {r,t)dVt 2 , 



(1) 



where dVl 2 is the line element for the unit 2-sphere. We 
use subscripts to denote partial derivatives. The Einstein 
field equations yield 



p 



m(r) 
R 



R 2 R.r'' 



(2) 
(3) 



the latter equation defining the density p(r, t) of the dust; 
m{r) is arbitrary (but subject to certain conditions in- 
troduced below) and will be referred to as the mass, al- 
though in fact m equals twice the Misner-Sharp mass. In 
particular, this locates the apparent horizon at R = m. 
Solving (||) gives 



R 3 {r,t) = -m(t c (r 



(4) 



where t c (r) is another arbitrary function which, in the 
process of collapse, corresponds to the time of arrival of 
each shell r to the singularity. 

We assume that the collapse proceeds from a regular 
initial state; i.e. at time t — (i) there arc no singu- 
larities (all curvature invariants are finite) and (ii) there 
are no trapped surfaces i.e. m(r) < R(r,0). We have the 
freedom of a coordinate rescaling in (Q) which we can use 
to set R(r, 0) = r. This gives 



, s 2 /r 3 
3 V m 



(5) 



Notice that this leaves just one arbitrary function, m(r). 
There is a curvature singularity called the shell-focussing 
singularity along t — t c (r), so the ranges of the coordi- 
nates r, t are < r < oo and < t < t c (r). The possi- 
bility that the shell-focussing singularity may be naked 
was first noted in ]T3] | , and the first mathematical results 
on the role of initial data were given in 0. We note 



This so-called shell-crossing singularity is gravitationally 
weak though what this means in terms of con- 

tinuation of the geometry is not yet known. Since we 
are primarily interested in the shell-focussing singularity 
we impose the condition that along each world-line r = 
constant, the shell-crossing singularity does not precede 
the shell-focussing singularity. That is, t c (r) < t sc (r) 
for all r > 0. This is equivalent to taking R r > for 
all r > 0, and yields rm >r < 3m, r > 0. We note that 
M0) = h. 

The initial data consists of just one function a(r) := 
p{r, 0) defined on an interval [0, b] where b is the initial 
radius of the collapsing dust sphere. The result of 0] 
confirmed in Q is as follows: 

Theorem 1 Let b > and n e C 3 [0, b]. Let 



m(r) 



x 2 u(x)dx 



defined on [0, b] satisfy the no-shell crossing condition 
rm. r < 3m on (0, b] . Then the marginally bound col- 
lapse of the dust sphere with initial radius b and initial 
density profile /x(r) results in a singularity from which 
there emanates a f.p. radial null geodesic if and only if 
one of the following conditions is satisfied. 

1. //(0) < 0. 

2. //(0) = and u"(0) < 0. 

3. //(0) = n"(0) = and 

u'"(0) 2 

< -3(26^/3 + 45) ~ -60.0222, 

where u := /x(0). 

In g], we found it convenient to make the following def- 
initions which will be used below. 



( \ 3 /MO . x 

m(r) =: r (— + mi), 

m tr (r) =: r 2 (p + m 2 ), 

m, rr (r) =: r(2u + m 3 ), 



mi(0) = 0; (7) 
(8) 

m 3 (0) = 0. (9) 



Notice that m 2 (r) = u(r) — u , and so m 2 G C 3 [0, b] and 
7712(0) = 0. These terms are related as follows: 
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Proposition 1 m\,mz € C 3 [0, b] and 



m'i(0) = -m' 2 (0), 
m'/(0) = im£(0), 
<(0) = \<(0), 



m' 3 (0) = 3m£(0); (10) 
m(f(0) = 4m£(0); (11) 
m' 3 "(0) = 5m' 2 "(0). (12) 



As we shall see below, the proofs of the results in ||] 
sometimes require a different approach in each of the fol- 
lowing cases: 

(i) m[(0)^0. 

(ii) mi(0) = 0,m'/(0) 7^ 0. 

(iii) mi(0) = 0,mi'(0) = 0,mi"(0) 7^ 0. 

III. TIME-LIKE GEODESICS 

In this section, we determine the circumstances under 
which radial and non-radial time-like geodesies may em- 
anate from the singularity. 

The governing equations for the geodesies with angular 
momentum L obtained from the Euler-Lagrange equa- 
tions are 



t + R r R rt r 2 + R,, f — 
R 6 



i? r f + R, rr r 2 + 2i? rt f i - 
-i 2 + (i? r /f + 



R :i 
I? 
R 2 



0, 



= 0, 



(13) 
(14) 
(15) 



where the over-dot represents differentiation with re- 
spect to an affine parameter s along the geodesies. We 
have e = +1 for time-like geodesies and e = for null 
geodesies. We are looking for the existence of a solution 
of (p^)-(p^|), which satisfies the following condition: 

Existence condition: 

There exists 8 > such that i is a non-negative, 
integrable function of the affine parameter s and r is an 
integrable function of s for s £ [0,5) and such that 



lim r(s) = 0, 



lim t(s) = to- 



Non-negativity of t implies that the geodesic is future- 
pointing. We shall now demonstrate the following result: 



Proving one half of Theorem 2 is quite straightforward. 
Considering the f.p. causal geodesies through a point p 
of this spherically symmetric space-time as curves in the 
r — t plane, we can compare slopes at p using ( |l5| ) to show 
that a f.p. causal geodesic 7 through p, which is not the 
unique outgoing radial null geodesic 7* through p, locally 
precedes 7*. That is if 7,7* are given by t = t 1 (r), 
t = i 7<i (r) respectively with i 7 (ro) = i 7 „ (ro) where r$ 
corresponds to p, then t lr (r) > i 7 (r) for r < r . Hence if 
7* avoids the singularity i.e. reaches r = at some time 
i* < t c (0), then so too must 7. A general result of this 
nature is given in p7| . So we have 

Proposition 2 // there are no f.p. radial null geodesies 
emanating from the singularity, then there are no f.p. 
causal geodesies emanating from the singularity. 

Proving the converse result is more difficult. However 
it turns out that the proof of Proposition 9 in 0] relating 
to non-radial null geodesies requires only minor modifi- 
cations in order to be applied in the present case. The 
idea of the proof is to identify a region Q bounded by a 
line r = 7'o > and by curves t = t\(r),t = t2{r) which 
intersect at the singularity and out of which time-like 
geodesies moving into the past may not pass. This con- 
finement is proven by carefully controlling the derivative 
of x/y, where x — t and y — R, r r, and the region may be 
constructed only if there are radial null geodesies ema- 
nating from the singularity. Using the geodesic equations 
(|l|)-(|l|), we obtain 



dfx 
ds\y 



L 2 

R 3 y 



G 1 




y 



(16) 



where G{r,t) = RR rt /R. r . The definition of f2 guaran- 
tees that G is positive therein. The key part of the proof 
for null geodesies involves showing that d(x/y)/ds > in 
n. But positivity in the null case (e = 0) clearly implies 
positivity in the time- like case (e = +1), and so the re- 
sult of Proposition 9 of ||] carries through. Thus we can 
state: 

Proposition 3 Let /i(r) be such that there exist f.p. ra- 
dial null geodesies emanating from the singularity. Then 
there exist f.p. radial and non-radial time-like geodesies 
emanating from the singularity. 

This completes the proof of Theorem 2 and, together with 
the results of (|], gives a complete account of the causal 
geodesies emanating from the shell-focussing singularity 
in marginally bound spherical dust collapse. 

In the next sections, we consider some implications of 
the existence of this array of geodesies emerging from the 
singularity. 



Theorem 2 Given regular initial data for marginally 
bound spherical dust collapse subject to a condition which 
rules out shell-crossing singularities, there are radial and 
non-radial timelike geodesies which emanate from the 
ensuing singularity if and only if there are radial null 
geodesies which emanate from the singularity. 



IV. STRUCTURE OF THE SINGULARITY I: 
CAUSAL STRUCTURE 

One more item is required to determine the causal 
structure of the singularity: the behaviour of ingoing ra- 
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dial null geodesies (IRNG's). This behaviour is the same 
for every case which gives rise to a naked singularity and 
is summarised as follows: 

Proposition 4 Let m be such that the initial data give 
rise to a naked singularity. Then there exists a unique in- 
going radial null geodesic which terminates at t = to,r = 
0. 

Proof: To see that such geodesies must exist, consider 
an outgoing radial null geodesic 7 which emerges from 
r = at some time t\ < t . Then the radial ingoing null 
geodesies originating on 7 reach r — either before, at 
or after t = to. Let Si be the set of values of r on 7 for 
which the first outcome holds, and let S2 be the set for 
which the last outcome holds. Then any IRNG through 
points of 7 for which r satisfies supjSi} < r < infjS^} 
must terminate at r = 0, t = to- 

In order to prove uniqueness, let e > and consider 
two IRNG's ti(r),t2(r) terminating at r = 0,t = t and 
satisfying ii(r) < t2(r) on (0, e). Note that the geodesies 
cannot intersect in t < t c (r). The equation governing 
IRNG's is 

— — -R — 1 Q-l/3 m - r *sc — t , 

dr~ ' r ~ m 2/3( tc _ t )i/3' 
and using t\{r) < fa{r) < to, we obtain 

J:te-*i)< 12 " 1/3 ^s (tc-*,)- 1 ' 8 ^-*!). 

If //(0) and /x"(0) are not both zero, then the coefficient 
of ti — 1\ on the right hand side is integrable on [0, e) and 
so £2(0) — ii(0) = gives t2 — <i = 0. In the case where 
/Lt'(O) = //'(0) = 0, fi"'(0) ^ the transformation u = r 3 , 

y = ~ m'" W ( i | 1 ) 3/ ' 2 (*0 — *) casts tne equation in a form 
such that exactly the same argument may be used. • 

The fact that the visible portion of the singularity may 
be identified with a point (0, to) on the boundary of a 
local coordinate chart means that this portion of the sin- 
gularity must be null. The fact that there is a unique 
ingoing radial null geodesic terminating at this portion 
of the singularity indicates that it has the form of an ingo- 
ing rather than an outgoing null hypersurface. Thus the 
local causal structure of the singularity is as indicated in 
Figure 1. The global structure of the space-time depends 
on the initial radius of the dust sphere, which will deter- 
mine whether the singularity is locally or globally naked. 
That is, if the initial radius of the dust sphere is suffi- 
ciently small, then future pointing radial null geodesies 
which emerge from the singularity will enter the vacuum 
region prior to the formation of the apparent horizon. 
Since the apparent horizon meets the event horizon at 
the boundary of the star, such geodesies must therefore 
lie in the portion of the Schwarzschild space-time out- 
side the event horizon, and extend to Sf + , giving a glob- 
ally naked singularity. There are also cases where all f.p. 
causal geodesies emerging from the singularity cross the 
apparent horizon and terminate at the future singularity. 
In those cases the singularity is locally naked. 



1 



<S : t = t c {r),r > 




FIG. 1: Local causal structure of the singularity in spheri- 
cal dust collapse. We identify three different portions of the 
singularity; Vo, the end-point of the time-like geodesic To 
(r — 0,t < to); A/0 (t — to,r = 0), the visible portion of the 
singularity and So (t = t c (r),r > 0), the censored portion of 
the singularity. iVo is the unique ingoing radial null geodesic 
terminating at To- 



V. STRUCTURE OF THE SINGULARITY II: 
SOME USEFUL LIMITS. 

In this section, we consider the limits of various quan- 
tities along different causal geodesies terminating at and 
emanating from the singularity. 



A. Causal geodesies running into the singularity 

We consider first the congruence of future-pointing 
time-like geodesies r = r > with tangent u = . This 
is the congruence formed by the world-lines of the dust 
particles and t is proper time along each member of the 
congruence. We study the approach to the singularity 
from the point of view of the behaviour of the character- 
istic lengths l a , a — 1,2,3 along the eigen-directions of 
the expansion tensor 9 a b = V a Ub, defined by 

' Q — a 
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where 9 a is the eigenvalue associated with the direction 
a. The geodesies in question have eigenvalues 



R,rt 
R r 



1 



{tc-ty 1 - {t sc -t)- 1 



and 



-ft,* 
R 



'-.(tc-t)- 1 . 



Then the characteristic lengths l a along r = satisfy (by 
virtue of t c (0) = t sc (0) = to) 

li=h=h = (to-t) 2/3 . 

Thus this singularity is a 'point' in the Thorne classifica- 
tion (see e.g. chapter one of ||). On the other hand, the 
characteristic lengths along r — tq > satisfy 



h = 



h = 



(t sc {r ) ~ t){t c (r ) ~ ty 1 ? 3 

(t.c(n>)-*c(n>))(tc(n>)-t)~ 1/3 , 
h = {t c {r Q )-tf'\ 



so these singularities (i.e. the different endpoints of these 
geodesies) are 'cigars' in this classification. Apart from 
the issue of visibility, this is the first indication of a sig- 
nificant distinction between the singularities Vq and Sq. 

The shear and expansion scalars of this congruence are 
given by 



f 



1\2 



((t sc -t)-' -(t c -t)-') 



o = ((tsc-tyi + itc-t)- 1 ), 

while the density and Weyl invariant are 
9 = \{tc-t)-\t sc ~t)-\ 

* 2 = ^(tc-ty^tsc-ty^tc-tsc). 

The limiting behaviour of ratios of these terms have 
significance in the study of singularities, in particular 
in the context of the Weyl curvature hypothesis and 
the study of isotropic singularities which has grown up 
around this. See |1 for the ideas which motivate this 
study, and JlSfl for a more recent review of the topic. The 
hypothesized scenario which these studies seek to corrob- 
orate or deny is that the initial cosmological singularity 
has low 'gravitational entropy' J8| and consequently sub- 
dominant Weyl curvature, while singularities arising from 
the end-state of collapse are highly entropic and so are 
accompanied by dominant Weyl curvature. 

We recall that isotropic singularities have been put for- 
ward as a model of the initial cosmological singularity, 
and the sub-dominance of the Weyl curvature is mani- 
fest in these models in that 



si r^abed 
yy2 WfccrfL- 



RefR £f 







(18) 



in the limit as the singularity is approached [91. Other 
limits may also be derived for the case of isotropic singu- 
larities M, for example 

,A a 2 n . A 

where A = G a bU a u h and u a is tangent to a time-like 
congruence. The limit is taken along this congruence, 
which is assumed to have a certain degree of regularity. 

Now, in the present case, along the central world-line 
To, we find the following limits: 

which reveals an interesting comparison with isotropic 
behaviour summarised above. In fact, the equality holds 
identically for the first three limits. The sign difference 
in the last term arises from the fact that we are studying 
collapse. 

Along the world- lines r = ro > 0, the corresponding 
limits are 



lim {W 

-»tc(r Q ) 



.4 

ff2' 



1 



-,6} = {+oo,0, -,-oc} 



Note that the shear is divergent in the approach to the 
singularity and that the first of these limits indicates the 
dominance of the Weyl curvature at Sq. The latter actu- 
ally holds for any causal curve running into So- To see 
this we can write 



1 1, 
W = -- 



tr 



6 t c -t 



Since t sc (r ) > t c (r ) for r > 0, this quantity diverges 
to +oo in the approach to the singularity t = t c (r), r > 
regardless of the direction of approach. 

Again, we conclude that there is a significant difference 
between the behaviour at Vo and at Sq. 

As far as ingoing radial null geodesies are concerned it 
was already demonstrated in section IV that there is a 
unique ingoing radial null geodesic which terminates at 
Vq. Repeating the arguments of the next section (from 
( |l9| ) onwards), we can shown that the Weyl-to-Ricci ratio 
W is finite in the limit as the singularity is approached 
along this geodesic. This ratio diverges in the approach 
to So along ingoing radial null geodesies. 



B. Causal curves emerging from the singularity 

As we have shown above, the singularity is naked if 
and only if there are time-like geodesies emerging the 
singularity. These results only show existence of such 
geodesies and give very little information about them, 
certainly not enough to carry out the analysis done in 
the previous section. However, there is one situation in 
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which this analysis is possible, namely self-similar col- 
lapse. This requires the use of a slightly different gauge 
from the one used above. The line-element is given by 
(U), but the solution for R is now []l9| 

R = r(l- Ax) 2 / 3 , 

where x — t/r and A is a constant parameter. Note that 
t = is no longer a regular initial data surface, but any 
t — ti < is such. The shell-focussing singularity occurs 
at x = A -1 . The singularity (r. t) = (0,0) is naked if 
and only if < A < A c ~ 0.64 p9[ . In the case where 
the singularity is naked, the surfaces x =constant are 
space-like for x < Xq and for x > Xi, and time-like for 
xq < x < Xi, where xq,xi are respectively the smallest 
and largest positive roots of x — i? r (x) =0. x = xq 
is the Cauchy horizon, i.e. is the first outgoing radial 
null geodesic emerging from the singularity. All outgoing 
radial null geodesies in the region xq < x < xi,r > 
originate at the singularity (0, 0). Thus for x c £ (xo, x±), 
the curves t = x c r, (#,</>) ^constant are time-like curves 
originating at the singularity. We analyse this congruence 
under the same headings as the previous section. 

Writing R, r = (1 - Xx)- 1 ? 3 ^ - Ax) = F(x), this con- 
gruence has unit tangent vector field 

u = a(x)(x- + -), 

where a(x) = (x 2 — F 2 )^ 1 ' 2 . The acceleration vector 
field is 

a(x) , d d N 
a = — - — - ( — — +x- 



dt dr 



where 



a(x) = -^X 2 a 2 F(l-Xx)- 4 / 3 . 

The eigenvalues of the expansion tensor are found to be 
all equal and are given by 

a(x) 



Any given member of this congruence has, from the form 
of u above, 



dr 

rf7 



a(x c ), 



and so the proper time r satisfies r = r / a{x c ). Thus the 
characteristic lengths all satisfy 



giving l a — t for a — 1,2,3. Again we see that a visible 
portion of the singularity consists of a singular region 
which are 'points' in Thome's classification. 



As the expansion tensor has three equal eigenvectors, 
the shear of this congruence vanishes identically. The 
Weyl-to-Ricci ratio is 



* 2 1 
P 3 



Ax)~ 



which is non-zero and finite in the limit as the singularity 
is approached along x = x c . ^ and p both individually 
diverge at the singularity. The expansion 9 = 3a(x)/r 
diverges. 

In the general (i.e. non-self similar case), we can an- 
alyze this last limit for the time-like geodesies emerging 
from the singularity whose existence has been proven in 
Proposition 3 above. The limit we wish to determine is 



Wo = 



\]/ 9 

hm|^| 



lim 



1 tsc tc 

S6 t P - 1 



where the limit is taken along the geodesic emerging from 
the singularity. The value of this limit may be obtained 
using rHopital's rule: 



W = lim 



1 1' - t 



r-»0 6 t' 



(19) 



where t' is the slope of the geodesic in the r ~ t plane. 
Referring to the proof of Proposition 9 in we recall 
that along such geodesies, the slope t'(r) of these curves 
in the r — t plane satisfies 

R,r < t' < {l + S)R. r , 

for some S > and sufficiently small. The significance 
of R r is that it gives the slope of the outgoing radial 
null geodesies and therefore the bounds above allow us 
to focus on the radial null geodesies. We know that such 
geodesies emerge into the region to < t < tn{r), where 
t = tff(r) = t c — |m is the apparent horizon. Using 
the latter bounds on t, we can obtain bounds on R r (see 
equation (p"7|)), i.e. on t' . These in turn may be integrated 
to improve the original bounds on t. Focussing on the 
case where m' x (0) ^ and carrying out this process twice 
yields bounds of the form 

cir 2/3 + m(r) < t' < c 2 r 2/3 + u 2 (r), 

where < c\ < c 2 and Ui = o(r 2 ' 3 ) for i = 1,2. We can 
then take the limit in ([HI) and find that 



1 1[ 
Wq = lim - 



t' 



06 t' - 1' 



In the case m[(0) — 0,to"(0) ^ 0, the limit gives 
Wo = 2/9. These limits are the same for every causal 
geodesic emerging from the singularity. Finally, for the 
case m'^O) = to"(0) = 0, m'"(0) ^ 0, the existence of the 
limit is proved in the same way as above. In this case, 
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the result depends on the value of the third derivative 
and on the geodesic in question. 

Our principal conclusion is thus that the visible por- 
tion Ao of the singularity does not exhibit the Weyl- 
dominance which is thought to be characteristic of singu- 
larities arising in collapse. Nor does it exhibit the Weyl- 
sub-dominance (18) characteristic of isotropic singulari- 
ties. In the case where the calculation may be carried 
out, the singularity is point-like. 

As an aside, we remark that we can find examples 
where the quantity W 2 increases or decreases with time 
depending on the geodesies under consideration. This is 
potentially of interest to debates on gravitational entropy 
(see e.g. §, M and 63) . 



C. Redshift 

In order to make contact with recent results given in 
H we study briefly the redshift along outgoing geodesies. 

Consider a photon which is emitted by a source with 
unit time-like tangent !t" s j and received by an observer 
with tangent u?-. . Let k a be the tangent to the photon's 
world-line. Then the redshift z relative to the emission 
and reception events P\ and Pi is given by 



[kaU a {s) } Pl 



[ k aU a [o) ]p 2 ' 

The redshift of photons emerging from a naked singular- 
ity has been calculated in [p2f and |3| for various trajec- 
tories. The source and receiver are taken to be co-moving 
with the fluid: take the source to be r = ro << 1 and 
the receiver to be r = r\ > r§. Thus the source hugs 
the central world-line To and the null singularity Ao and 
terminates at Sq near to the junction of Ao and Sq. Then 



(s) 



— °0' 



and 1 + z = tp 1 /tp 2 . Taking Pi to be 
the past endpoint of a geodesic emerging into the future 
from the singularity places the photon source on the sin- 
gularity and can be thought of as allowing ro | 0. In [ p2| , 
it was shown that the redshift thus measured is finite for 
singularities of classes (i) and (ii) of section II, but infi- 
nite for singularities of class (iii). In ||, the redshift is 
shown to be infinite for photons propagating along non- 
radial geodesies emerging from the singularity. This lat- 
ter result is immediate as we see from ([15]) : i diverges at 
the singularity R = for any non-radial causal geodesic 
which is not a radial null geodesic. This shows that the 
result - infinite redshift -also applies to particle de Broglie 
wavelengths (k a time- like). This is a possible (classical) 
indication that particles emerging from a naked singu- 
larity may not be associated with the transfer of infinite 
amounts of energy. Furthermore, if physical particles es- 
caping from the singularity always travel along geodesies 
with some angular momentum, then they will always be 
infinitely redshifted to an outside observer. We empha- 
size that these results are purely classical and it would 



be interesting to study them using e.g. a semi-classical 
approximation. 



VI. STRUCTURE OF THE SINGULARITY III: 
TOPOLOGY 



In the previous section, we classified the different por- 
tions of the singularity using Thome's point-cigar-barrel- 
pancake classification. This describes the behaviour of a 
co-moving fluid element carried by a time-like congru- 
ence meeting the singularity. However, this does not ad- 
dress another question which relates to the topology of 
the singularity and which we shall now describe. In the 
conformal diagram of Figure 1, each point in the two- 
dimensional representation of the non-singular part of 
the space-time manifold represents a 2-sphere: the point 
(r, t) of this diagram represents the 2-sphere with radius 
R = R(r, t) at time t. The question we address here is 
the following: does each point on the singular bound- 
ary Vo U Ao U So represent topologically a 2-sphere or a 
point? Christodoulou has shown that the future outgo- 
ing null (and hence censored) singular boundary which 
can arise in the collapse of a self-similar scalar field is fo- 
liated by points, while the space-like singular boundary 
that arises in another sector of the space of solutions is 
foliated by 2-spheres ||. 

Christodoulou's argument runs as follows. Any ra- 
dial null geodesic of the space-time actually represents 
a 2-sphere's worth of null geodesies. Consider a radial 
null geodesic 7 running into a future outgoing null sin- 
gular boundary Ai. This actually represents a family 
7*,£ G S 2 of radial null geodesies running into the sin- 
gularity. From this family, choose a fixed member 7^ 
and the antipodal member 7 s (N, S for north, south) . 
Now let {Pn' S } be causally increasing sequences of points 
along j N ' S respectively, chosen so that the singularity 
is reached in the limit n — ► 00 and for each n, the 
points p^' S are antipodal on the same 2-sphere S 2 . 
Christodoulou shows that the causal pasts J~(Pn' S ) of 
Pn' S satisfy 



lim J'{pZ) = lim J-(p b n ), 



and so the ideal boundary points lrnin^ooP^' , having 
the same causal past, are, by definition, identical. See 
Figure 2. 

The key ingredient required for this argument is to 
show that past-directed non-radial null geodesies through 
Pn orbit the axis R = sufficiently quickly (in the limit 
as n — ► 00) so that the causal past of p^ is contained in 
the causal past of ■ That is, the rapidly orbiting non- 
radial null geodesies have enough 'time' to get around 
from the north to the south pole of S 2 and so have a 
look at what is to be seen from there. See Figure 2. 

We argue here that, at least in the self-similar case, the 
topology of the naked portion of the singularity is not 
uniform, being point-like on one region and spherical on 
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X\ are outgoing radial null geodesies emerging from the 
singularity (t,r) = (0,0). We recall that x corresponds 
to the Cauchy horizon, i.e. the boundary of the future of 
a regular initial data surface t = ti < 0. 

ORNG's initially in the region x £ (xq,xi) remain 
in this region and satisfy x'(r) = (F — x)/r < 0. Thus 
along such a geodesic, x increases as r J, 0. But x is 
bounded above by X\, and so must have a finite limit 
as r { 0. This limit (lim- = lim^p) must solve x = F 
and so must equal x\. A similar argument applies for 
ORNG's in x e (x 1 ,X^ 1 ). Hence all ORNG's emerging 
from the singularity, except for the single ORNG x = xq , 
originate at (0,0) with x = x\. However as we will 
see, non-radial null geodesies may emerge with different 
values of x in this limit. We note that this indicates 
that considering only radial null geodesies does not give 
a complete picture of the singular boundary. 



FIG. 2: Behaviour of non-radial geodesies near the singu- 
larity. For convenience, we consider a future outgoing null 
singular boundary Mo. R = along Mo- The straight ar- 
rows represent ingoing radial null geodesies j N ' s with initial 
points Pq' S which are antipodal points on an initial 2-sphere 
So- These geodesies project to the same radial null geodesic 
in the Lorentzian 2-space dfi = 0. The curved arrows are 
past-directed non-radial null geodesies emerging from differ- 
ent points £ 7 , the different lengths of these curves rep- 
resenting the fact that angular velocity <j> of these geodesies 
increases as their initial point p 1 ^ approaches the singularity 
R = 0. In the case where (j> —* — oo as the singularity is 
approached, the geodesies reach all the way around the reg- 
ular center R = to the opposite side, and in the limit, we 
may consider that the non-radial geodesic emitted from p^ 
reaches p^. Hence as points of the causal boundary, p^, and 
pfj, are identified and Sec is a single point. In the case where 
<fi has a finite limit as the singularity is approached, Soo is a 
2-sphere. 



another. This requires a fairly thorough analysis of the 
behaviour of null geodesies emerging from the singularity. 

In the self-similar case, the line element may be 
written (see e.g. jH|) 

ds 2 = -dt 2 + F 2 (x)dr 2 + r 2 U 2 (x)dn 2 , (20) 

where x — t/r and 

U(x) = (1-Ax) 2 / 3 , 



A 



F(x)^(l-Xx)- 1 ^(l^^x). 



The shell- focussing singularity occurs at x = A -1 . 
The equations governing outgoing radial null geodesies 
(ORNG's) may be written as 

dx F — x 
dr r 

Then it can be shown that the singularity is naked for A 
in the range (0, A c ), in which case F — x has two positive 
roots xq < xi < A . The curves x = xq and x = 



We now consider the question raised above. Each 
of the radial null geodesies emerging from (0, 0) rep- 
resents a 2-sphere's worth of ORNG's emerging from 
the singularity. To address the issue of the topology of 
the singularity, we consider non-radial null geodesies 
(NRNG's) emerging from the singularity. 

The homothetic Killing vector field which generates 
the self-similarity of (|20|) is 



i-t— +r— 

dt dr 

If k a is tangent to a null geodesic, then g a b£, a k b is constant 
along the geodesic. This is a first integral of the null 
geodesic equations, which we may therefore write as 



i 2 + F 2 r 2 



L 2 



-t 



F 2 



r 2 U 2 
k 



(21) 
(22) 



where L is the conserved angular momentum and k is 
the constant generated by £. We can assume that k > 0, 
as a change of sign in k corresponds to a change in sign 
of the affine parameter. Our principal result which is 
of significance for the topology of the singularity is the 
following: 

Proposition 5 (i) Through every point of {(t,r) : xq < 
x < Xi}, there exists a non-radial null geodesic which 
emerges from the singularity (0,0). Along this geodesic, 
<t> is a non-integrable function of the affine parameter r , 
and (j> diverges to — oo in the limit as the singularity is 
approached. 

(ii) Through every point of {(t, r) : x\ < x < A -1 }, 
there exist non-radial null geodesies which emerge from 
the singularity (0,0). Along every such geodesic, <f> has a 
finite limit as the singularity is approached. 

The proof is given in the appendix. Proposition 5 
allows us to deduce the following results regarding the 
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topological nature of the singularity Ao , from the causal 
boundary perspective. 

We start by recalling that radial null geodesies emerge 
into the region xq < x < X\. Now, on any such geodesies, 
we may choose a point p N sufficiently close to the singu- 
larity such that there exists a non-radial null geodesic 
through p N with the property than this NRNG reaches 
the antipodal point of the 2-sphere (on which p N lies) 
before the radial geodesic reaches the singularity. The 
significance of this is that in the limit as the point p N 
approaches the singularity, this point and its antipodal 
partner p s must be considered to have the same causal 
future, and so are identified as points of the causal bound- 
ary. Therefore, the corresponding section of Ao is foliated 
by points. 

On the other hand, the results obtained for NRNG's 
in the region x\ < x < A -1 indicate that the NRNG's 
through a point p N do not have enough time to reach 
the antipodal point p s before the singularity is reached. 
Hence in the limit as the singularity is approached, the 
causal futures of p N and p s will not match up with one 
another, and hence these points cannot be identified as 
points of the causal boundary. In fact p N will not have 
the same limiting causal future as any other point on the 
same 2-sphere, and so the corresponding portion of Ao is 
foliated by 2-spheres. 

Finally, one further point of relevance may be deduced 
from the proof of Proposition 5: In every case considered, 
the future evolution of the NRNG's terminates at x = 
A -1 at some r > 0. The limiting behaviour of f and x in 
terms of x can be read off from the governing equations, 
and can be used to obtain the limiting behaviour of R = 
rU + rU'x. This can be integrated, and generically yields 
a result indicating integrability of 1/R 2 . Hence (f> has a 
finite limit along any NRNG approaching the singularity 
x = A -1 at some r > 0. Thus the singularity Sq is 
foliated by 2-spheres. 



VII. CONCLUSIONS AND COMMENTS 

We have extended the result of Q to all f.p. causal 
geodesies by proving that for the marginally bound spher- 
ical dust collapse there are radial and non-radial timclike 
geodesies which emanate from the ensuing singularity if 
and only if there are radial null geodesies which emanate 
from the singularity. It remains to be seen whether this 
result is true in general in spherical symmetry, at least 
for singularities at zero radius. The fact that the naked 
singularity always has time- like geodesies emerging there- 
from leads to interesting possibilities in terms of the spec- 
trum of particles that may be created by such a singu- 
larity: such particles may be massive as well as massless. 
This would seem to be vital if the information lost in the 
process of collapse is to re-emerge at a later stage in the 
form of a black hole or naked singularity explosion [^3| . 

We have studied the structure of the singularity from 
several points of view. We started by studying the 



behaviour of limits of various quantities along causal 
geodesies terminating at the singularity and emanating 
from the singularity. In this way we have classified dif- 
ferent sectors of the singularity as points and cigars. We 
have then considered the topological properties of the 
singularity. These studies were motivated by the fact 
that the existence of non-radial null geodesies emerging 
from the singularity suggests that it must have some non- 
pointlike structure. Naively, one might expect that the 
singularity Ao of Figure I must by cylindrical, i.e. fo- 
liated by 2-spheres, since points thereon emit a full 2- 
spheres worth of non-radial null geodesies. This is in 
contrast to a regular center of a spherically symmet- 
ric space-time, through which non-radial null geodesies 
cannot pass. However detailed calculations indicate that 
while the topology of So is that of a 2-sphere, the topol- 
ogy of Ao (in the sense of ideal boundary points) may be 
either that of a 2-sphere or of a point. 

The calculations in section V of various limits at the 
singularity seem to us to be of interest beyond merely 
classifying the different portions of the singularity. The 
question of what constitutes a genuine space-time singu- 
larity is not yet settled. This is of particular relevance 
for cosmic censorship hypothesis: a visible singularity 
which is mild in some suitably (and rigorously) defined 
way should not be considered a counter-example to the 
hypothesis. The following two definitions have proved to 
be useful. Both involve the question of predictability to 
the future of a singularity. 

We mention first that of Clarke [^4j. Put simply, a 
singularity is considered inessential if it does not present 
any obstruction to the evolution of test fields in space- 
time. The question of global hypcrbolicity of the space- 
time (cosmic censorship) is redirected to the question of 
generalized hyperbolicity of fields on the space-time. Ex- 
amples of space-times admitting singularities in the clas- 
sical sense of C 2 ~ geodesic incompleteness but satisfying 
a generalized hyperbolicity condition have been given in 
@, @ and {§. 
The second definition is that of an isotropic singularity 
. As mentioned above, this was introduced in the study 
of the initial cosmological singularity and hypotheses re- 
lated to this [|| . The connection with cosmic censorship 
is perhaps not immediately obvious, but becomes so when 
one acknowledges that the initial cosmological singularity 
is naked: every past-directed causal geodesic meets the 
singularity. The isotropic singularity program has pro- 
duced the following profound result: the big bang can be 
considered a regular initial data surface for the Einstein 
field equations. (Different matter models have been stud- 
ied; polytropic perfect fluids in j27j, the Einstein- Vlasov 
model in [g8| . These rely on results for singular differen- 
tial operators due to Newman-Claudel p9|).) Thus with 
isotropic singularities, we see that it is possible to have 
global hyperbolicity of the space-time to the future of a 
singularity. 

So with a view to addressing the question of whether or 
not the naked singularities studied here should be consid- 
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ered genuine counter-examples to cosmic censorship we 
could ask: do we have any evidence that the singularity 
is either (i) inessential in Clarke's sense or (ii) isotropic? 
We will not address the first of these here. Regarding the 
second, we contend that there is some evidence that, at 
least in the self-similar case, the singularity A/o U Vo has 
characteristics in common with isotropic singularities. 

As we have seen, the Weyl-sub-dominance property 
( |l8| ) - a characteristic of isotropic singularities - is not 
satisfied for the congruence of time-like curves emerg- 
ing from the singularity which we studied in Section 5.2. 
Nevertheless, we speculate that the fact that the ratio W 2 
is finite (but not zero) may turn out to be an indication of 
good behaviour of Mq regarded as an initial data surface. 
As evidence in favour of this, consider the following de- 
tails which relate again to the self-similar case. From the 
calculations in the appendix, we know that outgoing ra- 
dial null geodesies emerging from the singularity fall into 
two classes. The first class contains a single member, 
the ORNG x = xo, which emerges from Vq in Figure 1. 
The second class consists of all other ORNG's emerging 
from the singularity: these emerge from Ao • Along these 
geodesies, the limit x — > x\ is satisfied as the singularity 
is approached. Hence the congruence of time-like curves 
x = x c G (xq,x\) - call it the similarity congruence - 
considered in §V.B must be considered to emerge from 
Vo- This resembles the behaviour of the initial singular- 
ity in certain FLRW universes, namely those filled with 
a perfect fluid satisfying a equation of state p = a p. with 
a e (—1,-1/3]. In these space-times, the entire fluid 
congruence emerges from a single point of the singularity 
and a null portion of the singularity forms a portion of 
the causal future of this point. See [Q for details. This 
resemblance is borne out by writing the line-element in 
coordinates adapted to the similarity congruence. Intro- 
duce a time function T = rH{x) chosen so that the sim- 
ilarity congruence is orthogonal to T =constant. Then 
the line-element becomes 



ds 2 



{x 2 - F 2 ) 



( x 2 - F 2 ) 2 ^- + F 2 dx 2 + U 2 dQ'' 



metrics mentioned above, and is similar to the decompo- 
sition whose existence defines isotropic singularities (the 
difference is that t ^ -oo at the singularity currently 
under discussion). 

So while there is still no direct evidence that this naked 
singularity does not break the future predictability of the 
space-time, the question appears to be worthy of future 
investigation. 
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APPENDIX A: PROOF OF PROPOSITION 5 

(i) Given x c S (xo,Xi), it is straightforward to show 
that there are NRNG's confined to this time-like hyper- 
surface which have the desired properties. Substituting 
x =constant into (21) and (E2h yields 



= ±^(x 2 - F 2 )- 1 ' 2 = k(a 



F 2 )-\ 



This applies for r > and x € (xq, x\), wherein x 2 - 
0. The function H satisfies 



F 2 > 



which is valid since x 2 > F 2 in (xq,X\). Choosing k, L so 
that these two match, which gives a 1-parameter solution, 
gives the geodesic in question. A sign choice is required 
in order that the geodesic is future pointing t > 0. Note 
that rf — c, which is positive for the future-pointing 
choice and so r 2 ~ 2cr as r J, 0, choosing the origin of 
the affine parameter t to coincide with the singularity. 
Then 



w_ 
IT 



p2 ' 



and therefore is finite on (x$,x\). The line-clement may 
be written as ds 2 — T 2 ds 2 , where 

ds 2 = -a 2 (x)^- + b 2 (x)dx 2 + c 2 (x)d0 2 , 
1 1 

with a, 6, c finite and non-zero on (xo,xi). T is the ex- 
ponential of proper time along the similarity congruence. 
This decomposition of the metric into an overall scaling 
factor (which vanishes at the singularity) times a non- 
singular metric (when written in terms of r = In T) is ex- 
actly the same decomposition which exists for the FLRW 



r 2c ' 

and so <f> diverges to — oo as the singularity is approached, 
(ii) To prove this part of the proposition requires the 
study of all non-radial null geodesies in the region 
(x±, X" 1 )- The equations ( pi] ) and (B2T) yield a quadratic 
equation for r, which has solutions given by 



1 



X " 2 i L , 2 



F 2 

These give 



-\-k±-(k 2 + —(F 

: L p y U 2 



x 2 )) 



l/2l 



T 2 

— IF 2 

U 2 1 



x 2 )) 



1/2 



(Al) 



(A2) 
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The solutions with the lower sign require the change 
t — * — t in order to be explicitly future-pointing. Imple- 
menting this, the equations are (any two of) 



1 ^[ fc +|( fc2 + §2( F2 -- 2 )) l ]>°( A3 ) 



X 



F 2 - 
1 

r 2p 



(k 2 + ^(F 2 - x 2 )) 1 / 2 >0 



F 2 k 
t = — r > 0. 

x rx 



(A4) 
(A5) 



Recall that k > 0. Positivity of x immediately gives 

t < xiT (A6) 

for any point (t, r) preceding a chosen initial point (ij, rj). 
Positivity of f and k yields 



dx 
dr 



< 



F 2 



< 2(F[-l)(x-x 1 ) + M(x-x 1 ) 2 



There latter inequality arises from using Taylor's theo- 
rem on F, which is analytic at x±. 2M > is a bound 
on the closed interval [a;i,a;,] of the (continuous) second 
derivative of (F 2 — x 2 )/x. F[ = F'(xi), and this number 



is always bigger than 1. Integrating over 
< r < ri yields 

t > xir + ( 



ci 



1 — c 2 r- 



-)r 



A+l 



where 



(A7) 



where C1.2 a re p ositive con stants and A = 2(F[ — 1) > 0. 
Combining ( |A6| ) and (A7), we see that these geodesies 
all originate at [t = 0,r = 0). 

The rate at which r shrinks to zero can be obtained as 
follows. We have 



2k 



rr 



F 2 



while 



x ~ — r 2 . 

X\ 

These may be integrated to give 

r ~ fcrV^i'+l) = kr a 
where a < 1/2 and so 

L 







dv 



converges to a finite value in the limit as r \ 0. 

Next we consider th e s olutio ns c orresponding to the 
choice of upper sign in (Al) and ( A2). We note that t > 
and x > are automatic consequences of this choice. On 
the other hand, f > if and only if H{x) > 0, where 



H 



L 2 
U 2 



fc 2 



We note that H'{x) > 0. Then for any initial value H(tq) 
of H , the increase of x brings about the increase of H to 
positive values. Then r becomes and remains positive, 
and so the geodesic must run into the singularity x — 
A -1 . The past evolution requires more careful attention. 
We deal with the case k > first. 

We have rr — rr(x), and by using a Taylor expansion 
centered on xi, we find 



rf{x) 



Hjxi) 
2k 



+ 0{x-xi). 



If in the past evolution we reach t\ < tq whereat H{t{) < 
and x(ti) > xi, then H < for all r < n. Hence r 
increases into the past, and so the geodesic inevitably 
extends back to x = x\ at some r > 0, thus avoiding 
the singularity. This situation is obtained if and only if 
H(x 1 )<0. 

So it remains to consider the cases H{x\) > 0. If 
H{x\) > 0, then the argument above shows that rr > 
for all x > x\ . 

We can obtain the bound rx' (r) > (F 2 — x 2 )/x. Inte- 
grating over [r, r,-] (0 < r < r;) yields 



^dx>\n\ r A. 



F 2 



Thus if the geodesic approaches r = in the past, then 
the integral on the left hand side here must diverge. This 
can only occur if we also have x | x\ as the singularity 
is approached. Let this occur at r = t s . But then 



X\ 



= lim - 

dt 

= hm — 

t-^o dr 

F 2 

= lim ( — 

l^H x 



xi 



xrr 

k 2 L 2 k 2 , 
2 -(772 " Zl)" 1 , 



xi U 2 



which yields a contradiction. Therefore these geodesies 
must extend back to x = Xi before reaching r = 0. 
This argument applies when H(xi) > and requires only 
slight modification to be applied to the case H(x\) = 0. 
Note that there is a finite (rather than infinitesmal) in- 
terval of affine parameter between an initial point on the 
geodesic and the point corresponding to entry into the 
region x <x\. 

Finally, we need to consider the NRNG's with k = 
0. Note that in this case, the choice of sign in (Al) 



and (A2) is only of a time-orientation and so without 
loss of generality we only treat the choice of upper sign. 
The argument leading up to (A7) may be repeated to 
show that in this case, there exist geodesies emerging 
from (0, 0), and that x — > x\ along these geodesies as the 
singularity is approached. It remains to determine the 
rate at which r approaches zero for these geodesies. 
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From the governing equations, we have x = G(x)r 2 , 
where 

We can also write 

(r 2 y =2r 2 x(F 2 -x 2 )- 1 x. 
These yield the second order equation 

x = J(x)x 2 , 

where 

,G 2x 



we obtain 



J(x) = ( 



G F 2 



r) 



2 l F{ - 1 



)(x — Xi) 1 + 0(1) as x — > Xi 



This equation may be integrated to obtain the 
r— dependence of X cLS X I x\. Then using r 2 = H/x, 



C 2 T 2/(F[ + 1) 



Recalling that F[ > 1, we see that the power here is less 
than one. Hence 



L 



L 



R 2 U 2 { Xl )r 2 



is integrable in the limit as the singularity is approached, 
giving a finite limit for the value of <f>. 

To summarise: NRNG's in x > Xi either exit this re- 
gion at some finite radius, taking a finite amount of time 
to do so, or emerge from the naked singularity (0, 0). In 
the latter case, <f> has a finite limit in the approach to the 
singularity. This concludes the proof of the proposition. 
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